Antiferromagnetism and Superconductivity in a model of Quasi ID Organic 
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I apply a two-step renormalization group method to the study of the competition between antifer- 
romagnetism (AFM) and superconductivity in an anisotropic 2D Hubbard modeL I show that this 
simple model captures the essentials of the ground-state phases of the quasi ID organic conductors. 
As found experimentally, the ground-state phase diagram is mostly AFM. The AFM is localized 
in the strong-coupling limit where the electrons are confined in the chains. It is an SDW in the 
weak-coupling limit where interchain hopping is present. There is a tiny region in the weak-coupling 
regime where transverse two-particle hopping is dominant over magnetism. 



Introduction. The intriguing discovery of supercon- 
ductivity (SC) lying next to antiferromagnetism (AFM) 
in the phase diagram of the charge transfer Bechgaard 
salts (see FiglT]) remains one of the greatest issues of 



condensed matter physics The proximity of AFM 



and SC turned out to be a generic feature not only of 
the Bechgaard salt series {TMTSF)2X, but also of the 
Fabre salts (TMTTF)2X and layered 2D organic and 
cuprate superconductors. In the quasi ID organic mate- 
rials, AFM occupies a large region of the phase diagram 
and is believed to be central to the emergence of SC. It is 
believed that the understanding of AFM is a prerequisite 
to that of SC. 

The nature of pairing in these compounds is still un- 
clear. Recent experiments have yielded conflicting re- 
sults. A NMR Knight shift experiment by Lee et al. fi^ 
found that the symmetry of the Cooper pairs is triplet in 
{TMTSF)2{PF)c. No shift was found in the magnetic 
susceptibility at the transition for measurement made un- 
der a magnetic field of about 1.4 Tesla. A subsequent 
Knight shift experimencnt performed at lower fields re- 
veals a decrease in the spin susceptibility. This result is 
consistent with singlet pairing. [16] The authors of this 
later experiment suggested a possible singlet-triplet pair- 
ing crossover as function of the magnetic field as a reso- 
lution of these conflicting results. 

In the face of these experimental uncertainties, a theo- 
retical input onto the behavior of simple models of these 
compounds is of crucial importance. A theoretical anal- 
ysis of AFM in the quasi ID organic materials was pro- 
posed by Bourbonnais, Caron, and coworkers [l[. This 
description was essentially based on a perturbative renor- 
malization group (RG) applied on the g-ology model. 
They found that the AFM phase has two regions. On 
the left side of the AFM phase, the ID chains are Mott 
insulators, the charge gap Ap induced by coulomb inter- 
actions is such that Ap » t±^ , where t±^ is the transverse 
hopping parameter. Hence, the single particle transverse 
hopping is irrelevant. In this region, the electrons are 
necessarily confined in the chains. t±^ can nevertheless 
generate an interchain exchange Jl by virtual interchain 



hopping. Using the RG method , it can be shown that 
this process leads to a transverse effective Hamiltonian 
i/_L = / dxY,i Ja.S{x)iS{x)i+i, "^itl^ '^J- ~ *i/^p- In 
this region, the electrons are necessarily confined in the 
chains due to the irrelevance oit±. As pressure increases, 
the electrons progressively delocalize in the transverse 
direction. In the right region of the AFM phase, the 
magnetism is itinerant and arises from the nesting of the 
Fermi surface Cr^iik) — — e;.j.(A: -t- Q), where the indices 
(r, I) stand for the right and the left parts of the Fermi 
surface respectively, and Q = (2A:f, tt) is the nesting vec- 
tor. The nesting leads to the divergence of the suscepti- 



/(£,-(fc))-/(e,(fc+Q)) 



However, 



bility x(Q, w = 0) oc ^ ,^(fc)„,,(fc+Q) 
this description presents a difhculty. At a temperature 
above the SDW phase, the system is not a Fermi liquid, 
it is rather a Luttinger liquid. There are no quasi parti- 
cles and the FL description does not apply. A perturba- 
tive RG applied for coupled LLs can only yield the most 
divergent susceptibility, but it cannot reach the ordered 
phase. Further increasing pressure destroys the magnetic 
order and leads to superconductivity as illustrated on the 
right part of the phase diagram. A recent RG analysis 
has suggested that superconductivity can emerge from 
the frustration induced by hopping between next-nearest 
neighbors chains [2^. Pairing emerges through short- 
range AFM analogous to the Kohn-Luttinger effect in- 
duced by Friedel oscillations [Hj. 

Model and Methods. The Hubbard model is also of- 
ten used for the organic conductors. It differs from the 
g-ology model in that it includes non-linear dispersion 
and the scattering processes are not restricted to a nar- 
row region around the Fermi points. While it enjoys 
an exact solution and can be studied with the density- 
matrix renormalization group (DMRG) or the Monte 
Carlo methods in ID, it remains a difficult challenge, even 
numerically, when ij_ is turned on. In this letter, I will 
show that the essential aspects of the ground-state phase 
diagram can be obtained from the single-band extended 
Hubbard model using the two-step DMRG method [6, ^, 
to which I supplement with a Wilson RG analysis of the 
low energy spectrum. I will consider the following model 
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FIG. 1: Sketch of generic phase diagram of the quasi ID 
organic conductors: LL (Luttinger Liquid), FL (Fermi Liq- 
uid), MI (Mott Insulator), CO (Charge Ordering), SP (Spin- 
Peierls), AFM (Antiferromagnet), SC (Superconductor). 



at quarter filling, the nominal density of organic conduc- 
tors, which has standard notations: 

H = -ijj ^{cl i^^Ci+i^i^o- + h.c.) + U ^n,,ijni,i,4 

X(c^^^^^Cij+i,o- + h.c.) + VI X n.i^i,<Tni,i-\-i^a- (1) 

i,l,<7 i.l,a 

The two-step RG method starts by using DMRG to 
compute the low energy eigenvalues and eigenvectors of 
a single chain of lenght L. Then Hamiltonian H is pro- 
jected onto the tensor product of the eigenvectors of the 
disconnected chains. This leads to the effective one- 
dimensional Hamiltonian 

^ ~ ^ Ha^i - t± ^(cJo-Ci+i.o- + h.c.) 

I 1,(7 

+Kl ^ "-i,i,o-ni,i+l,cr, (2) 
i,l.{T 

where i?o,i is diagonal, its elements are the eigenvalues of 
the single chain, the operators ci^^r are composite matrices 
containing the renormalized Ci^i^a for i = 1,L. Since H 
is 1£>, it may be studied using DMRG again. 

Additional insight in the behavior of H may be gained 
by using the Wilson RG instead of DMRG in the second 
step. The advantage of the Wilson RG lies in the fact 
that the low energy spectrum can be obtained. However 
the Wilson method, directly applied to i?, is not accurate 
because all the terms in the transverse direction are of 
the same order. I use the same trick used by Wilson for 
the Kondo problem [2^. H is defined as the limit of H\ 
when A — > 1. H/^ is given by 



I 

where 

+V1 ^ ni,;,crni,i-(-i,o-- (4) 

It is to be noted that this is not the usual Wilson's 
momentum space discretization. The justification of this 
scheme rests on the fact that the essential physics remains 
unchanged. When A > 1, the terms corresponding to 
I > 1 act as a perturbation on the term with ^ = 1. For 
A not too large, I expect H\ to essentially have the same 
behavior as H. But if A ^ 1, i_L/A will be too small 
with respect to the finite size energy separation and the 
chains will be disconnected. It is to be remarked that 
this approach may also be useful if it is embedded as a 
cluster solver in a chain-dynamical mean-field approach 

a. 

Results. The pressure variation will be mimicked by 
varying the Coulomb parameters U and V, while I keep 
1 11 , t± and V± constant in most simulations. In the regime 
of strong U and V, the isolated chains are Mott insula- 
tors, there is a large charge gap Ap while the spin degrees 
of freedom are gapless. At I^ ^ 2 for any value of C/, there 
is a insulator-metal transition [2^. The transverse cor- 
relation functions yield information on which type of or- 
der will dominate. Since the parallel correlations have a 
power law decay in absence of a gap, the most dominant 
transverse correlation in a given channel will automati- 
cally lead to long-range order in that channel, even if it 
is not dominant in ID. Before performing the analysis of 
the transverse correlations, it is somewhat instructive to 
look at the low energy spectrum provided by the Wilson 
method. This will allow to make a qualitative compari- 
son with the evolution under pressure with the prediction 
of the perturbative RG. 

The low energy-excitation spectrum of H was obtained 
from the Wilson RG for A « 1.2 for a lattice size 16 x 6, 
keeping 100 states block. The lowest 1000 excited states 
are shown in FigO They are drastically different as pres- 
sure is varied. In the confined regime for U ~ 6 V = 2 
FigUl all the lowest states have the charge A„ = 0. This 
is consistent with the fact that since t± ^ Ap, the low 
energy behavior of H is roughly identical to that of an 
Heisenberg model, only spin excitations are allowed. This 
is typically the regime of Fabre salts at ambient pressure 
where a large char ge g ap is observed in optical conductiv- 
ity measurements 21]. By contrast, in the SDW regime 
for U = A V = 0.85, spins excitation are still the lowest 
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FIG. 2: Charge An of lowest 1000 excitations after 6 RG iter- 
ations, the lattice size is 16 x 6, for A = 1.2: in the localized 
AFM [/ = 6, V = 2(left), the SDW C/ = 4, V = 0.85 (center), 
and superconductor U = 2, V = (right) 



but excitations with A = ±1 now appear above them. Ex- 
citations with A = ±2 are also observed at higher energy. 
It is to be noted that because of numerical errors, exci- 
tations with A„ = 1 and A„ = — 1 which should normally 
have the same energy are shifted. In the regime with 
important superconductive correlations U — 2 V — 0, 
excitations with A = ±2 now appear closer to the ground 
state. This is consistent with the increase of supercon- 
ductive correlations as we found below. 

I now analyse the evolution of the transverse corre- 
lations when U and V are varied using the two-step 
DMRG for the lattice size x Ly = 16 x 17. I keep 
msi = 256 states during the first step and a maximum 
of ms2 = 128 states during the second step. For this 
value of ms2, AE/t± « 5 which means that we are at 
the limit of the two-step method as discussed in Refd. 
Starting from the left of the AFM phase where U and 
V are expected to be strong, because of the presence 
of a large Ap, the carriers are confined in the chains, 
even when tj_ <C Ap is turned on. The carrier con- 
finement was observed by Vescoli et. al in optical re- 
flectivity measurements 22]. When the oscillating elec- 
tric field was oriented in the transverse direction, no 
plasma mode was observed. This carrier confinement is 
seen in the behavior of the transverse Green's function 
G{y) = (cL/2,L/2+ycl,/2,L/2+i)- ^(y), shown in FiglSta) 
for [/ = 6 and V — 2, decays very fast. G{y) « for 



y > 3- This was expected given that tj_/Ap « 0.1. But 
as predicted by the RG [H, although irrelevant, t± can 
nevertheless generate the motion of transverse spin de- 
grees of freedom and lead to magnetic order. This is seen 
in the transverse spin-spin correlation function C'{y) — 
|(SL/2,L/2+ySi/2x/2+i) which is showu in FiglSl^b). I 
find that despite the irrelevance oft±, G{y) has its largest 
amplitude in the strong coupling limit. As expected I find 
that the transverse singlet (triplet) superconductive cor- 
relations SSiy) = 2(Ai/2,L/2+yA[/2,L/2+i) iST{y) = 

2(0L/2,L/2+y0i/2,L/2+l))' '^^'^^^ = ^ (Ci,iT C,+i,;x - 

Ci,iiCi+ij{) {Qij = ^(cij|Q+i,^ arc neg- 
ligible in this limit as seen in Fig[3l^c)(d). 

Moving toward the right of the phase diagram by in- 
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FIG. 3: Transverse interchain correlations as function of dis- 
tance y: Green's function G{y) (a), spin-spin correlation C{y) 
(b), singlet superconductive SS{y) (c), triplet superconduc- 
tive ST{y) (d) for [/ = 6, V = 2 (circles), f/ = 4, V = 0.85 
(squares), U — S, V = (diamonds), U — 2, V = (trian- 
gles). 



creasing pressure or equivalently decreasing U and V, 
the carriers are expected to progressively deconfine. In 

re- 



this regime, the reflectivity measurements of Refl22 
ported the observation a transverse plasma mode. The 
carrier deconfinement is expected for t±/Ap w 0.5. For 
U — 4 and V — 0.85, the numerical simulation yields 
tj_/Ap « 0.3 In Figl3]Ja),(b),(c),(d), it can be seen that 
for [/ — A and V = 0.85 G{y) now has a non-zero am- 
plitude, C{y) is still significant, while SS{y) and ST{y) 
are still very small. This suggests that the system is in 
the SDW phase. The perturbative RG, coming from high 
temperatures, shows that there is a ID to 2D crossover 
at Tx w t±/Tr. At Tx the ID RG equations cease to be 
valid. FL arguments are used to describe the onset of 
the SDW order. The experimental observations are how- 
ever that both the non-FL and FL characters seem to be 
present depending on the quantity measured A fu- 
ture application of the TS-DMRG at finite temperature 
could shed light on this ineteresting crossover regime. 

If U and V are further reduced, C(y) now decays faster 
despite the fact that the amplitude of G{y) is larger. C{y) 
now vanishes for y > 3 as seen in Figl3Ub) for U = 3, 
U — 2 and = in both cases. This implies the ab- 
sence of long-range magnetic order in this regime. At the 
same time, SS{y) sharply increases suggesting the onset 
of superconductivity in agreement with the phase dia- 



4 



gram. The values of SS{y) at long distances are however 
within our margin of error. The presence of V± is crucial 
to the enhancement of pairing correlation. In the absence 
of Vj^, if all other parameters are unchanged, the domi- 
nant correlations are AFM. ST{y) is still negligeable thus 
suggesting the absence of triplet pairing in the extended 
Hubbard model. Experimental results on the symmetry 
of the pairs are still controversial. Knight shift experi- 
ments from different groups have predicted singlet [l6l | 
and triplet pairings as discussed in the introduction. 
This study shows a tendency towards singlet-pairing only 
in the extended Hubbard model. Singlet pairing, though 
interchain, was also predicted by the perturbative RG. 
However, in the RG, the pairs could be formed by carriers 
lying on neighboring chains and AFM was suppressed by 
hopping to next-nearest neighbor chain. It was suggested 
in the RG study that the interchain pairing could be 
triplet in presence of strong enough next-nearest neigh- 
bor hopping and nearest neighbor Coulomb interaction 
both in the transverse direction. The present two-step 
results do not however settle this issue. Many small ef- 
fects including longer range hopping and Coulomb inter- 
actions were not included in this simple model. Triplet 
superconductivity could emerge from these terms. 

Conclusion. To summarize, I have shown that the sin- 
gle band extended Hubbard model displays the ground- 
state phases of the quasi ID organic conductors: (i) lo- 
calized magnetism in the strong-coupling regime, (ii) de- 
localized SDW magnetism in the intermediate coupling 
regime, and (iii) possible superconductivity of singlet 
type in the weak coupling regime. I have not discussed 
the spin-Peierls phase which rests at the extreme left of 
the phase diagram. In this regime, the electron-phonon 
coupling is dominant over the effective transverse ex- 
change Jl, hence a pure ID study of a spin model coupled 
to phonons such as that of Ref. 24] , which shows a spin 
gap opening, captures this part of the phase diagram. 
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